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For a long time, one of the main limitations of cesium atomic fountains has been the cold collision 
frequency shift Q jl^llOtlTlLte^ . By using a method based on a transfer of population by adiabatic 
passage (AP) H71 Il9| allowing to prepare cold atomic samples with a well defined ratio of 
atomic density as well as atom number we have a better measurement of this effect. By improving 
the method we found out some unexpected properties of cold collisions, confirming that 133 Cs is a 
rich atom to study cold collisions. Those results lead to fine comparisons with cold collision theory 
and may constraint some parameters of cesium. 
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Frequency standards are used in many industrial 
(Global Positioning System, navigation) and scientific 
pursues (fundamental physic tests [5|, Il5l |l6( ) . Primary 
standards are today based on cesium and are contribut- 
ing to the realization of the second. The most accurate 
are using cold atoms to allow long interaction time. In 
cold atomic gases collisions play an important role by 
shifting the clock frequency. Thus, it is crucial to control 
and to study accurately cold collisions for clock uncer- 
tainty budget. Nevertheless, there is a tradeoff between 
the clock stability and its accuracy. In this paper we will 
describe improvement of the a technic based on fast AP 
that allows to reconcile those two goals. For high den- 
sity atomic sample, such as in Bose-Einstein condensates 
(BEC) 1], manifestations of collision such as molecular 
Feshbach resonances are observed. We found out, 

in our fountain, such phenomenon at very low magnetic 
field and density. Theoretical tools will be remind and 
a numerical simulation of some of our data will be pre- 
sented. 



THE DOUBLE FOUNTAIN: RECENT 
IMPROVEMENTS AND PERFORMANCES. 



Atom loading improvement. 

At Paris observatory, there is a set of three atomic 
fountains, using cesium. In particular, our one is a dou- 
ble fountain 0^0, able to work alternatively with 
rubidium and cesium. In the last year, on the cesium 



part, we improved the loading rate by supplementing a 
transverse collimation of the atomic beam in addition to 
the present chirp cooling (see FIG. It consists in a 
2D optical molasses with two laser beams, recycling the 
light in a "zigzag" way j^,!^. This allows us to cool 10 9 
atoms in 200 ms in an optical molasses. The 2D molasses 
increase the loading rate by a factor of ~ 10. Regular pa- 
rameters of the double fountain are a magnetic field in 
the interrogation zone of 2 mG, a loading duration of ~ 
200 ms that allows to load 10 9 atoms within a fountain 
cycle of ~ 1.2 s. The launch velocity is ~ 4 m.s -1 that 
imply a launch height of ~ 0.8 m. 



Frequency stability. 

A sapphire cryogenic oscillator (SCO) 13] from the 
University of Western Australia is now used as a filter for 
one of our H-maser. The SCO is an experiment by itself 
|2f)l] , but is very reliable and can be used routinely to drive 
the fountain. The SCO is weakly phase-locked on the H- 
maser with a time constant of about an hour. Therefore, 
the resulting signal reflects the good short term stability 
of the SCO while in long time scales it reproduce the 
characteristics of the H-maser. The SCO delivers a signal 
of about 12 GHz. 

At BNM-SYRTE most of the microwave synthesizer 
use a 100 MHz signal to generate the 9.2 GHz necessary 
to interrogate the atoms. Therefore, a down-converter 
has been developed to transform the SCO signal into a 
100 MHz one. An evaluation has shown that the foun- 
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FIG. 1: Schematic of the double fountain both base on Rb 
and 133 Cs, so that there is two pairs of cavity finely tuned on 
hyperfine atomic frequency of corresponding atoms and two 
ovens. On the cesium part, a transversal collimation of the 
atomic beam has been improved in addition of the present 
chirp cooling. This consists in a 2D recycling light, pair of 
laser beam, used in a "zigzag" way. This allows to load around 
10 9 atoms in 200 ms in an optical molasses. 



tain was limited by the phase noise of the microwave 
setup. Henceforth, a new microwave synthesizer able 
to down-convert directly the 12 GHz to 9.2 GHz has 
been built in order to minimize the phase noise. It 
achieves a frequency stability of ~ 3 10~ 15 at Is. In 
order to take benefits of its excellent short term stability 
(~ 5 10~ 16 at Is) a large number of trapped atoms is 
required |2l|. With ~ 10 7 detected atoms we were able 
to reach 1.6 10~ 14 at 1 s, the best frequency stability for 
a primary standard to date. With this performances a 
resolution of 10 -16 is achievable in 6 hours. 

FIG.[21shows the Allan deviation of the fountain driven 
by the SCO. This is a differential measurement [23j, this 
means that a full atom number configuration (■) is al- 
ternated with a half atom number configuration (•) each 
60 fountain cycles. The (■) represents the fractional fre- 
quency deviation of the corrections applied to the interro- 
gating microwave generated with the SCO signal to stay 
at resonance with the atoms. This is performed with the 
maximum atom number {N at ) which allows the 1.6 10~ 14 
at Is stability. On short term (10 s to 100 s) the allan 
variance is dominated by the Cs fountain noise. After 
200 s, the drift of the SCO is observed. The Allan devia- 
tion of the frequency difference between low (N at /2) and 
high(A^ at ) density (A) averages as white frequency noise 
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FIG. 2: Shows the Allan deviation of the fountain driven by 
the SCO. This is a differential measurement of the collisional 
shift, this means that a full atom number configuration (■) is 
alternated with a half atom number configuration (•) each 60 
fountain cycles. The (■) represents the fractional frequency 
deviation of the corrections applied to the interrogating mi- 
crowave generated with the SCO signal to stay at resonance 
with the atoms. This is performed with the maximum atom 
number which allows the 1.6 10~ 14 at 1 s stability. The Al- 
lan deviation of the difference (A) averages well for all times 
and has also a good stability. This leads to a good statistical 
evaluation of the mean difference in a very reasonable time. 



for all times and reaches ~ 2 10" 16 for a 20000 s measure- 
ment. This leads to a good statistical evaluation of the 
mean difference in very reasonable time. The SCO drift 
is efficiently rejected and atom number dependent effects 
(i.e. cavity pulling Q and collisional shift) are measured 
with a statistical uncertainty of 10~ 16 in one day. This 
behavior indicates that determination of the collisional 
frequency shift is only limited by the frequency stability 
of the fountain. 



ACCURATE CONTROL OF COLD COLLISION 
FREQUENCY SHIFT. 

In previous work |l6l Il9j| , we demonstrated the possi- 
bility to control and evaluate the cold collision shift at 
the percent level by using the adiabatic passage method. 
Cold collisions represent the main systematic effect that 
shifts the frequency and restricts the accuracy in cesium 
fountains. At that time, our collisional shift for full den- 
sity, was around 10 -14 . Therefore, the corresponding 
uncertainty was near 10 -16 . 



Adiabatic passage performances. 

To achieve a stability around ~ 2 10 -14 at 1 s it is 
necessary to load a high number of atoms. For example 
the collisional shift corresponding to N at of the (■) curve 
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FIG. 3: Allan deviation of the ratio of atom numbers in the 
|4; 0} obtained by alternating each 60 fountain cycle a full 
atom number configuration over half atom number configura- 
tion on the |4; 0} state as a function of cycle number. The Al- 
lan deviation of the ratio averages as white noise and reaches 
~ 10 -3 resolution after 1 hour. It leads to evaluate AP at the 
10" 3 level. 



investigated in a future work. 

Similarly, we also check routinely the same ratio on the 
F = 3 level. Here, we find a larger deviation: 
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The departure of 3 1CP 3 from the 1/2 ratio is due to 
residual populations in the |3;771f) states. A fraction of 
the atoms is coming from the push beam (see FIG. ^ 
i.e. a laser beam that finalize to the selection process) 
which is de-pumping atoms from the |4; ttif) states to 
the |3;mp) states. The remanning fraction seem to be 
launched in the F = 3 state. Attempts to remove these 
atoms with extra repumper beams where partially un- 
successful. At this point, there is nothing much to do to 
prevent it but to look at the contribution of those spuri- 
ous atoms on the clock shift. 



Towards an adiabatic passage at the 10" 3 level. 



implies a frequency shift of about ~ 1.2 10 -13 . Thus, an 
AP at the percent level is no longer enough to stay at the 
10~ 16 level of uncertainty on this effect. To perform AP 
the amplitude of the microwave into the selection cavity 
has to be modulated and its frequency has to be swept. 
It is critical that the synthesizer used to chirp the selec- 
tion microwave stops the frequency sweep at resonance 
for the half AP. Our model (Marconi 2030) was not ac- 
curate enough and failed to be reproducible. It has been 
replaced by a DDS more easily frequency controllable. 
According to 0, 0] , one observable for AP evaluation 
is the ratio of atom number at full and half density. Each 
measurement is daily analyzed using this criterion to ver- 
ify that the AP method is properly working. 

FIG. |3 shows the Allan deviation of the ratio of atom 
numbers in the |4; 0} obtained by alternating each 60 
fountain cycle a full atom number configuration over half 
atom number configuration as a function of cycle num- 
ber. The Allan deviation of the ratio averages as white 
noise and reaches ~ 10 -3 resolution after 1 hour. It leads 
to evaluate AP at the 10 -3 level such as: 
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Where njvj? (n\^°l) is the density of the |4;0) atoms 
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for full atom number (half atom number) and A 
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(-Au.o\) is the detected atom number on the |4;0) for 
full atom number (half atom number). 

The deviation measured at the 10 -3 level could be in- 
duced by a non linearity in the detection. This will be 



An accurate control of the collisional shift is possible 
only if the contribution of the 0.3 % atoms populating 
states different from the clock states is correctly evalu- 
ated. To measure this, we intentionally populated the 
atomic cloud using AP with |3;tof ^ 0} states in addi- 
tion of the |3;0) clock state. An additional microwave 
synthesizer tuned on a |4; ^ 0} — > |3; m-p ^ 0) transi- 
tion is combined in the selection cavity with the regular 
selection synthesizer. A magnetic field is applied during 
the selection by AP in order to prevent parasitic excita- 
tions between zeeman sub-levels. To select, for example 
rnp = 1 state, one has to detune from resonance the ex- 
tra synthesizer of the corresponding first order Zeeman 
shift value. 

The measurement method alternate three different 
fountain configurations. The first one is a full AP selec- 
tion of the 3; 0) state, the second is a half AP selection 
of the |3;0) state. The first configurations give the col- 
lisional shift due to the clock states. The third one is a 
full AP selection of the |3; 0) state plus full AP selection 
of one of the |3;tof 0) states. The last configuration 
combined with the two others lead to a measurement the 
collisional shift due to the extra population. The respec- 
tive frequency shifts are: 
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Where noo («m F m F ) represents the atom density of the 
clock states (jhf states), and Kqq (K mFTnF ) is the cold 
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FIG. 4: displays J?om F as a function of the m F value for 
regular clock conditions. The contribution to the cold col- 
lision shift of the |tof| = 1; 2 are negative and equal for the 
same absolute values. At most the contribution of spurious 
|3;mp 7^ 0) is 1/3 of the contribution of clock states. The 
contribution of the Impl = 3 is different. 



collisional shift coefficient resulting of the extrapolated 
to zero density for the clock states (mp states). 

Under regular operating conditions, we evaluate the 
contribution to the clock shift due to the collisional shift 
of the extra |3; mp ^ 0) atoms per detected atom (v mF m F 
in Hz. (detected atoms) -1 ) normalized to the collisional 
shift of the clock states |3;0) atoms per detected atoms 
(foo) as a function of the mp value. This ratio (i?om P ) 
can be expressed as: 
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is the detected atoms number for 



the clock states (mF state) and can be expressed as 
iVoo = ^4|3 ; o) + ^4|4;0)- A verification that voo is field 
independent has been performed so that it could be used 
as a reference. It is important to ensure that the den- 
sity is the same whatever the could is made of. Thus a 
verification that the spatial distribution of the different 
Zeeman sub-levels is homogenous into the cloud has been 
performed. By selecting one m-p at the time and pulsing 
the Ramsey cavity with different durations, we access at 
different cloud slice. It shows that the peak density are 
equal to less than 7 %. 

FIG. 0] displays Rom F as a function of the tof value for 
regular conditions (B ~ 2 mG) . The contribution to the 
cold collision shift of the |mp| = 1; 2 are negative and 
equal for the same absolute values. At most the contri- 
bution of spurious |3; mp ^ 0) is 1/3 of the contribution 
of clock states. Thus, the collisional shift is controlled 
at the 10 -3 level. Surprisingly the contribution of the 
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FIG. 5: Ratio 7?oi as a function of the magnetic field. It 
represents a Feshbach resonance for [t7Zk | = 1. 



| Tof | = 3 is different, which was unpredicted at such low 
magnetic field. This observation is a clear indication of 
the presence of Feshbach resonances. 



FESHBACH RESONANCES. 
Experimental results. 

We thus looked for the ratio (i?om F ), for each absolute 
value of the mF as a function of magnetic field. We found 
out three Feshbach resonances (FIGs. El and 0) deeply 
field dependent. It gives a direct access to the amplitude 
and the width of each resonance. The field is very well 
controlled so that there is no way that error bars could 
mix up the shapes. For some values of the magnetic 
field the frequency shift due to |3;mp) =/= population 
can be as large as the standard clock shift. Let's remind 
that our regular magnetic field is 2 mG. One has to be 
very cautious by choosing a value of magnetic field in an 
experiment using cold 133 Cs. The |mp| = 2 resonance 
is obviously a multiple resonance, but the |tof| = 1 and 
specially |mp| = 3 look to be more simple. Assuming 
that | tof | = 3 resonance is a simple resonance we try to 
understand the line shape with a simple model. 



Theoretical approach: model for S-matrix. 

In order to describe the measured resonances, we need 
to account for the inelastic collision processes. Feshbach 
resonances in cold atomic systems are usually considered 
only in the context of elastic scattering, with one open 
channel only [Til Hif . In such a situation, only an elastic 
resonance width r e is involved, and the incoming chan- 
nel for the atoms is simultaneously the exit channel. Here 
we consider a situation with more than one open chan- 
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FIG. 6: Ratio i?02 as a function of the magnetic field. It rep- 
resents a Feshbach resonance for \m-p\ — 2. This is obviously 
a multiple resonance. 
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FIG. 7: Ratio i?03 as a function of the magnetic field. It repre- 
sents a Feshbach resonance for \m^\ = 3. This resonance has 
a more simple shape. Assuming that it is a simple resonance 
we gone to simulate the scattering process. 



nel. Atoms which exit in a different channel than the 
incoming channel will gain additional kinetic energy, and 
are considered as lost. This decay process can be medi- 
ated via the molecular Feshbach state, and gives rise to 
an additional energy scale, the inelastic resonance width 

The S-matrix for such an inelastic resonance can be 
written as 



S{k) = S bg (k) f 1- 



*r e (fc) 



E-u + Ur e (k) + Ur t (k) 



(i) 



with E = h 2 k 2 /m the relative kinetic energy between 
the colliding atoms, v = Afi(B — Bo) the detuning, and 
S bg (k) = exp(—2ikdb g ) the background (or direct) part 
of the S-matrix. Here Bq is the magnetic field value of 
resonance, and A/z the difference between the magnetic 



moments of the two particles in the incoming channel and 
the molecular Feshbach state. The background scattering 
length a6 9 = al g + ia bg contains not only a real part a£ g , 
but also an imaginary part a\ g since there is also a direct 
contribution to the inelastic decay. 

The energy widths are in principle a function of the 
relative wavenumber k, since they depend on the overlap 
between the molecular wavefunction corresponding to the 
Feshbach state, and the wavefunction in the correspond- 
ing open channel. In case of the elastic energy width, 
the energy dependence of the incoming channel gives rise 
to linear dependence in wavenumber, i.e. T e = c e k with 
c e a constant [l^. For the inelastic energy width, how- 
ever, the energy dependence of the corresponding open 
channels can be safely neglected over the energy range of 
interest. Therefore r, can be taken energy-independent. 



Monte Carlo model. 

The frequency shift of the clock transition due to the 
\3;mp ^ 0) state population is given by Eq. J5J): 
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where a and (3 refer to the |3;0) and 4; 0) clock states 
respectively. p 77 is the atomic density of the |3;tof 0) 
state (denoted as 7). To be consistent with our initial 
intent to analyze the |tof| = 3 data as a single Fes- 
hbach resonance, we will assume that a resonance oc- 
curs for only one of the two entrance channels (017 and 
/?7). The other entrance channel contribute only with a 
constant background part, independent on the magnetic 
field. This factor can be formally accounted for in the 
background term of the channel where the resonance oc- 
curs. Therefore, we are assuming that S' Q7 = 1 (resp. 
Sp-y = 1) while Sf3j(k) (resp. S a ^(k)) follows Eq. (Q. 

Eq. (|2J) clearly shows that the clock shift depends on 
the wavenumber k. The total clock shift is obtained by 
suitably averaging Eq. J3J) over the duration of the inter- 
rogation, taking into account the atomic cloud space and 
velocity distribution. This is done using a Monte Carlo 
simulation where the atomic trajectories are randomly 
sampled according to the space and velocity distributions 
as measured in the experiment (with CCD imaging and 
analysis of time of flight signals). The simulation ac- 
counts for the distribution of collisional energy, for the 
time varying (decreasing) atomic density, for the non- 
uniform response of the final transition probability to a 
time varying perturbation (sensitivity function), for the 
truncation of the atomic cloud crossing the microwave 
resonator. 

FIGs. 00 and ED represent examples of calculated line 
shape of the Feshbach resonance for various parameters. 
For all cases, we take a,b g — 0, B = and A/i = 
Also, we assume that the resonance occurs on the f3j 
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FIG. 8: For Ti = 0. For large values of C e , the shape and 
the ("natural") width of the resonance is determined by C e . 
The resonance has a symmetrical dispersive shape. For small 
values of C e , the shape and the width of the resonance is de- 
termined by the collision energy distribution. The line shape 
shows a strong asymmetry. The calculation clearly indicates 
that the top of the sharp feature corresponds to the magnetic 
field that meets the resonant condition. 
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FIG. 9: illustrates the effect of parameter Ti. C e is kept 
constant. The resonance shape is computed for different Ti 
amplitudes. The comparison with the measured line shape 
clearly indicates that the observed resonance are not strongly 
affected by the inelastic processes. The calculation shows that 
the sign of the clock shift is inverted when the resonance occur 
on the ay channel. 



channel. The initial distribution is gaussian with a = 
3.5 mm. The velocity distribution is oc (1+v 2 /v^)~ b with 
v c = 10 mm.s -1 and b = 2.1. The vertical axis represents 
the fractional frequency shift due to the 7 (|3; mp)) state 
for 10 8 atoms launched in this state. This corresponds 
to a typical effective atomic density of 2.2 10 7 cm~ 3 . 
In FIG. El T. t = 0. C e is equal to 0.2 x E rec /k rec , 
2 x E rec /k rec , 20 x E rec /k rec and 200 x E rec /k rec where 
E rec = h 2 k 2 ec /2m and k rec — 2-k/X, A = 852 nm. For 
large values of C e , the shape and the ("natural" ) width of 
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FIG. 10: The |?tof| = 3 data are in good agreement with 
a resonance dominated by the collision energy. The Monte 
Carlo model fits well the point and indicates the resonant 
field. 



the resonance is determined by C e . The resonance has a 
symmetrical dispersive shape. For small values of C e , the 
shape and the width of the resonance is determined by 
the collision energy distribution. The line shape shows a 
strong asymmetry. The calculation clearly indicates that 
the top of the sharp feature corresponds to the magnetic 
field that meets the resonant condition (vanishing differ- 
ence between the molecular bound state energy and the 
dissociation threshold). Similarly, FIG. E3 illustrates the 
effect of parameter IV C e is kept constant and equal to 
0.2 x E rec /k rec . The resonance shape is computed for r, 
equal to 0.02 x E rec , 0.2 x E rec , 2 x E rec and 20 x E rec . 
The comparison with the measured line shape clearly in- 
dicates that the observed resonance are not strongly af- 
fected by the inelastic processes. Finally, the calculation 
(as well as direct consideration of Eq. (JJJ shows that the 
sign of the clock shift is inverted when the resonance oc- 
cur on the ay channel. 

The measured |tof| = 3 resonance is fitted (see 
FIG. fTTTf> using the above model leaving B , A/i, C e , Ti, 
dbg as free parameters. The model also includes the pos- 
sibility of a slow variation of the background scattering 
length with the magnetic field through the addition of a 
linear and quadratic term as a function of the magnetic 
field. We checked that the following analysis is not sig- 
nificantly affected by this additional term. Finally, the 
model includes a scale factor for the atomic density. Sim- 
ilarly, the inclusion of this scale factor has no influence 
on the following analysis. 

The simple consideration of the sign of the data in- 
dicates that the resonance occurs in the (3y channel in- 
volving the upper clock state |4;0) and the |3;3) state. 
FIG. EH shows the optimized fit to the data. The qual- 
itative agreement with the model is quite satisfactory. 
For A/z, C e and Ti, the analysis does not lead to precise 
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value, due to the sensitivity of these parameters to the 
atomic cloud space and velocity distribution. We find 
A^i ~ 1.5jUb, C e ~ 0.2E rec /k rec and 1^ < 0.2 x E rec . 
Conversely, Bq is tightly constrained and well decorre- 
lated from the other parameters. Although an uncer- 
tainty on Bq cannot be easily extracted from our (non- 
linear) fitting procedure, we can safely state that Bq is 
constrained to within 1 mG (to be understood as a la 
error bar). We find Bq = 5 ± 1 mG, this corresponds to 
500 ± 100 nK, expressed in temperature scale. This is, to 
our knowledge, the lower molecular bound state energy 
ever involved in Feshbach resonances (i.e. around 10 kHz 
from the continuum). 

CONCLUSION. 

To sum- up, our primary goal is achieved, the colli- 
sional shift is controlled at the 10~ 3 of its value. A 
calibration of Zeeman sub-states contribution to the 
clock shift as a function of the field has been performed 
and teaches us magnetic field values to avoid. Under 
regular clock conditions an upper limit for their contri- 
bution shows an upper limit compatible with AP at the 
10 -3 level. This makes reachable the goal of primary 
standards at the 10~ 16 uncertainty level. Feshbach 
resonances have been observed for the first time at very 
low magnetic field and with a very good resolution. A 
Monte Carlo simulation has been performed and could 
fit properly some of experimental data. This constrains 
some parameters of the theory of collisions. 
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